Suppression of Quantum Fluctuations by the Dzyaloshinsky-Moriya Interaction in a 

Kagome Antiferromagnet 
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We report an electron spin resonance investigation of the geometrically frustrated spin-1/2 kagome 
antiferromagnet vesignieite, BaCu3V208(OH)2. Analysis of the line widths and line shifts indicates 
the dominance of in-plane Dzyaloshinsky-Moriya anisotropy that is argued to suppress strongly 
quantum spin fluctuations and thus to enforce long-range ordering. We also evidence an enhanced 
spin-phonon contribution that might originate from a lattice instability and discuss the origin of a 
low-temperature mismatch between intrinsic and bulk susceptibility in terms of local inhomogeneity. 
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The two-dimensional spin-1/2 quantum kagome anti- 
ferromagnet (QKA) has been extensively studied because 
of its geometrically frustrated magnetism [l|, |2|- The ex- 
act nature of its ground state (GS) remains mysterious 
even for the simplest nearest-neighbor isotropic Heisen- 
berg exchange (J), as several competingstates have been 
suggested - gapped [3|, |j| and gapless [5|, |6| spin liquids, 
as well as valence-bond solids [7|, |8| . Due to the appar- 
ent criticality of this model, further weak perturbing in- 
teractions give additional flavor to the subject. Indeed, 
the out-of-plane Dzyaloshinsky-Moriya (DM) magnetic 
anisotropy Dz [9J , that is present in any real QKA as the 
bonds lack inversion symmetry, is theoretically predicted 
to create a quantum critical point (QCP) at D'^/J ~ 0.1 
[10|. This separates a moment-free phase {D^ < D'i,) 
from a Neel ordered phase {D^ > D%) [l3-[i3- The min- 
eral herbertsmithite, 7-ZnCu3(OH)6Cl2, believed so far 
to be the best realization of the QKA [13|, appears to 
be a spin liquid |l4l . Il5| , in line with the value of its DM 
anisotropy, D^/J = 0.06(2) d, [13, [II- Its location in 
the region close to the QCP is likely to be related to ob- 
served field- [l9| and pressure-induced [20| freezing. On 
the other hand, the Neel ordered side of the phase dia- 
gram remains experimentally unexplored. 

The mineral vesignieite, BaCu3V208(OH)2, which 

crystallizes in the monoclinic space group C2/m [2]J , has 
been recently highlighted as a new realization of the QKA 
|22| . The minute 0.07% bond-length difference due to 
two inequivalent Cu^+ sites (Fi g, [p makes the triangles 
very close to being equilateral |23| . Indeed, there have 
been suggestions that the actual structure has equilateral 
symmetry |24| . though this may not yet be conclusive 
|25| . The magnetism of vesignieite is dominated by the 
nearest-neighbor antiferromagnetic interaction J = 53 K 
[22| that leads to a maximum in local susceptibility Xi 
at the temperature T ~ 0.5 J [26|, detected by nuclear 



magnetic resonance (NMR) . In marked contrast to her- 
bertsmithite [23|, [2J, |26|, [271 , vesignieite shows a magnetic 
transition to a g = Neel state at Tp{ = 9 K. During this 
spin freezing transition, an additional out-of-plane spin 
component creates a ZFC/FC bifurcation [27|. Based 
on the width of electron-spin-resonance (ESR) spectra 
|28| vesignieite has been suggested to possess large DM 
anisotropy D^ and thus to be in the ordered region of 
the phase diagram [23, l2J, [26|, [27|. Since vesignieite ap- 
pears to be the first clear case of a long-range ordered 
QKA and no proper attempt to identify and quantify its 
magnetic anisotropy has been reported, a detailed study 
is essential. 




FIG. 1: Two basic triangles with the two inequivalent cop- 
per sites Cu(l) and Cu(2) in the ab crystallographic plane 
of vesignieite. (a) The double-headed arrows connect apical 
0(1) sites in each CuOe octahedron and deflne the principal 
axis of the g-tensor on each Cu^^ site, (b) The Dzyaloshinsky- 
Moriya (DM) vector pattern of out-of-plane Dz and in-plane 
Dp components. In both triangles the spins are counted in 
the same circular direction leading to uniform D^. (c) Two 
principal directions, A and E, of the symmetric anisotropic- 
exchange (AE) tensor on each bond; A is canted by 6^0 = 45° 
out of the kagome plane, while E lies in the plane. 



In this Letter, we clarify the driving force of mag- 
netic ordering in vesignieite by determining its dominant 
magnetic anisotropy. Employing the local-probe ESR 
technique we show that the in-plane component of the 
dominant DM anisotropy, Dp, exceeds D^, in contrast 
to herbertsmithite. We propose that such a DM vector 
crucially suppresses quantum fluctuations and thus criti- 
cally affects the GS of this material. Additionally, we as- 
sess the importance of a symmetric anisotropic exchange 
(AE) that has recently been suggested as an important 
spin-Hamiltonian component of herbertsmithite [29|, ISOf . 

ESR has proven extremely efficient for determining 
magnetic anisotropy in spin systems, either by detecting 
collective excitations [3l|, l32|, or through the modeling 
of shifts [33|, |34| and line widths [5|, |35| of a paramag- 
netic resonance. Both, the shifts and widths are non-zero 
only when the anisotropy is finite |36| . and thus allow 
its direct quantification, as well as distinction of different 
forms of the anisotropy. Our experiments were conducted 
at 328.8 GHz on a custom-made homodyne spectrome- 
ter working in transmission mode at the NHMFL, Talla- 
hassee, USA. The sample was hydrothermally annealed 
powder similar to that used in previous studies |23l.l26|. 

In Fig. m^a) we show derivative ESR spectra typical 
of those recorded in the T-range between 3 and 300 K. 
The spectra have similar width as in herbertsmithite [5| , 
suggesting that substantial magnetic anisotropy is also 
present in vesignieite. In addition to this broad compo- 
nent, we observe a narrow component with the principal 
q-factors values of 2.05 and 2.25, typical of Cu^"*" ions 
[36|. We attribute this narrow component to a minor 
impurity phase since its intensity at 300 K amounts to 
only 0.3% of the broad-component intensity and exhibits 
a Curie-like T-dependence. The impurity signal is thus 
much too small to explain the substantial low-T increase 
of the bulk susceptibility Xb [Fig.IDJb)]. The ESR inten- 
sity of the broad component, XesRi convincingly follows 
Xb [Fig- E^b)], and not the nonmonotonic intrinsic Xi- 
This has an important implication for the hitherto un- 
known origin of the low-T increase of Xb ^^, |23| . Indeed, 
the observation of a single broad-component ESR line, 
rather than two distinct lines, reveals that the spins con- 
tributing to Xb are necessarily exchange coupled with the 
intrinsic Cu^^ spins. Bond disorder due to oxygen va- 
cancies or some other non-stoichiometry effect then pro- 
vides a credible explanation for the mismatch between Xb 
and Xi- Such disorder also explains the inhomogeneous 
broadening of the ^^V NMR lines far above T^ |26| . 

The experimental ESR spectra can be fitted [Fig.IJJa)] 
to a powder-averaged line shape based on a field 
distribution originating from the ji-factor anisotropy 
g(^0) = (g^cos^B + gpSin^Oy/^ that is convoluted with 
a Lorentzian function with the line width AB(9) — 
(AB^cos^e + AB^sin^Oy/^. Here z and p denote the 
orientations perpendicular and within the kagome plane, 
respectively, and 6 the polar angle between the applied 
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FIG. 2: (a) The ESR spectra (symbols) of vesignieite mea- 
sured at 328.8 GHz with fits (lines) further explained in the 
text. Arrows point to impurity lines with corresponding g- 
factors. (b) Comparison of the ESR intensity xesr (symbols) 
with bulk magnetic susceptibility Xfc (solid line) measured in 
5 T and the intrinsic susceptibility Xi (dashed line) obtained 
from ^^V NMR |26l |. The horizontal arrow indicates the in- 
tensity of the impurity ESR signal at 3 K. 



magnetic field and the normal to the plane. The pre- 
sumed independence of both g and AB on the azimuthal 
angle within the plane stems from the near threefold rota- 
tional symmetry of the spin lattice. Our approach elabo- 
rates on the previous ESR report that employed a simpler 
analysis yielding only an average AB and g [28] . We find 
9z > 9p- As the local anisotropy axis, set by the direction 
of the the shortened Cu-O(l) bond [Fig.[IJa)], makes the 
angle of 26.7° < 7r/4 with the kagome plane, the princi- 
pal g-iactor value gu along the apical direction is smaller 
than the value g± in the perpendicular direction. This 
confirms that the Cu^"*" orbital state involves occupation 
of the d322_r2 |22| rather than the more common d^^^-y^ 
orbital that would lead to g\\ > g± |36| . 

The ESR line width in vesignieite exhibits a minimum 
at Tmin = 40 K and increases linearly with T at least 
up to 300 K, i.e., T/J ~ 6 [Fig. [3ta)]. This is in sharp 
contrast to herbertsmithite where it was found constant 
for T > J [5]. In general, dying out of spin correla- 
tions above the characteristic exchange T causes a van- 
ishing contribution to the ESR line width. We therefore 
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FIG. 3: (a) T-dependence of the ESR line width (open sym- 
bols) at 328.8 GHz for the perpendicular (z) and the in-plane 
(p) directions of the magnetic field. The linearly increasing 
contribution (lines) is subtracted from the data to obtain the 
line width intrinsic to the kagome spin system (full symbols), 
(b) The AE model leads to a T- variation of the g-factor (lines) 
which is not detected in the experimental data (symbols). 



attribute the observed behavior to an additional line- 
broadening mechanism that arises from a spin-phonon 
coupling. A direct phonon process is expected to yield a 
linearly increasing ESR relaxation due to phonon mod- 
ulation of magnetic anisotropy [37|; e.g., as was found 
in SrCu2 (603)2 |35| where a structural phase transition 
was observed. We propose that in vesignieite the linear 
increase might be related to a structural instability as- 
sociated with the energetic proximity of the monoclinic 
crystal structure to the higher-symmetry rhombohedral 
{R3m) structure, and the perfectly undistorted kagome 
lattice [24| . Subtraction of the linearly increasing contri- 
bution gives the ESR line width intrinsic to the kagome 
spin system. Its increase with decreasing T below Tj^m 
is similar to that observed in herbertsmithite and can be 
attributed to the building-up of spin correlations, which 
are also responsible for the maximum in Xi at 25 K |26| . 
In order to determine the magnetic anisotropy, we there- 
fore make use of the 40 K spectrum, which corresponds 
well to the paramagnetic limit, as both the spin-phonon 
and the spin-correlation induced broadenings are small 
at this temperature. 
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FIG. 4: Reduced x^ of fitting the 40 K ESR spectrum with 
(a) the DM and (b) the AE model. The former yields opti- 
mal parameters \Dp\/J = 0.19, \Dz\/J = 0.07 and the latter 
two inequivalent solutions, A/J = ±0.15, E/J = =p0.13 and 
A/ J = ±0.04, E/J = =p0.21. Center of each panel: Compar- 
ison of the best fit (line) and experimental data (symbols). 



We model the ESR line-width anisotropy AB{9) by 
employing the well-established Kubo-Tomita moment 
approach [l|,|3|; valid in the paramagnetic limit. Both, 
the DM magnetic anisotropy H'y^j^ = J^fn) '^ij ' ^* ^ ^j 
[Fig[TJb)] and the traceless symmetric anisotropic 
exchange [FiglUc)], written in the local basis as Hae = 
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yield the line width of the general form 



AB{9) = V2tt 



kh 



{a + b cos^ 



2g{9)fiBJy c + d cos^e 



(1) 



where fc^ is the Boltzman constant, /i^ the Bohr mag- 
neton and the constants a, 6, c, d are relate to the 
anisotropy constants of each model (see Ref. |40|) . 

Fitting the experimental spectrum to the powder- 
averaged line shape with the line width given by Eq. ([6]) 
provides fits of equal quality for both models. These are 
displayed in Fig. [5] together with x^ maps spanned over 
the parameter space. For the DM model, we find the 
solution \Dp\/J = 0.19, |D^|/J = 0.07, while the AE 
model yields two inequivalent solutions, A/J = ±0.15, 
E/J = =f0.13 and A/J = ±0.04, E/J = ±0.21. The rel- 
ative sizes of the DM and AE anisotropy are very similar 



with respect to the dominant isotropic exchange J . Since 
the DM interaction results from a first order correction of 
J in the spin-orbit couphng, which is for the Cu^+ ions 
a ~^10% perturbation on J [36|, while the AE interac- 
tion is a second-order correction, the DM interaction is 
generally considered dominant. However, caution is nec- 
essary as Z?p, the largest anisotropy in the DM model in 
vesignieite, is reducible on the kagome lattice 10|. This 
is because it possesses a hidden symmetry [41| and can 
be transformed into an effective term of the order D'i/ J 
by applying a nonuniform spin rotation [42|. Therefore, 
we provide a second criterion that is based on the ESR 
line shift, which allows us to distinguish between the two 
anisotropy models in vesignieite. 

When the anisotropy is smaller than the isotropic ex- 
change, as is expected to be the case here, Nagata's the- 
ory [43| of the ESR line shift can be applied. Accord- 
ingly, the shift of the g-factor from its infinite-T value 
g°° is given by the first moment 43|, |4j|, g — g°° = 
([S*", [S+ ,n']]) /2^ibBq (S'^) where [] denotes a commu- 
tator, ( ) canonical averaging and S" = J2i ^i' the a- 
component of the total spin operator. It is important to 
stress that in this first-order calculation (in H') the DM 
interaction leads to zero shift |4J|, while the shift due to 
the AE interaction scales with the susceptibility x in the 
paramagnetic regime |43l . |45| , 
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Here A^^ denotes the component of the AE tensor along 
the applied field and A'^'^, A^^ in two perpendicular 
directions, 7V^ is the Avogadro number and \xq the vac- 
uum permeability. Quantifying this expression in vesig- 
nieite with the above-determined AE parameters, we ar- 
rive at the scalings {gz—gT)/x = ±22.3 (emu/mol Cu)~^ 
and [gp — g^)/x = Tll-7 (emu/mol Cu)~^ for the two 
relevant directions. These values significantly overesti- 
mate the measured ESR shifts that are found constant 
for T > 50 K within much smaller error bars [Fig. [3]Jb)]. 
The g^ data limit the AE anisotropy to at least 4-times 
smaller values. Therefore, the ESR line width, which 
scales with the square of the anisotropy, is entirely deter- 
mined by the dominant DM interaction. 

The experimental ESR shift becomes sizable below 
25 K [Fig. [3Kb)], which was ascribed previously to short- 
range ordering effects [28|. However, due to the traceless 
nature of the AE interaction it is generally expected that 
even in the short-range correlated regime the ESR line 
shift will exhibit both positive and ne gati ve shifts for 
different directions of the applied field [43 



45j . just as 



occurs at higher temperatures. This is not the case here, 
therefore, we attribute the solely positive low-T (/-shifts 
to local inhomogeneity present in vesignieite. Namely, in 
inhomogeneous systems (e.g., impure systems and sys- 
tems with several inequivalent sites) spatially varying lo- 
cal fields exclusively lead to positive g-shifts |46l . |47| that 



stem from the q^ Fourier components of the inho- 
mogeneous field |47| . This corroborates our proposition, 
based on the scaling of xesr with Xb, that there is sig- 
nificant bond disorder in vesignieite. 

Although, in both herbertsmithite [5| and in vesig- 
nieite the ESR spectra can be accounted for by the DM 
magnetic anisotropy, we point out an essential differ- 
ence. While the out-of-plane DM component D^ is dom- 
inant in herbertsmithite, it is the in-plane component Dp 
that dominates in vesignieite. We additionally note that 
the Kubo-Tomita approach generally underestimates the 
strength of the reducible Dp component with respect to 
the irreducible Dz component. 

Then, it is important to inspect the effect of both DM 
components on the GS of the QKA. Although in the clas- 
sical limit both components immediately lead to mag- 
netic ordering, their effect is rather different [48|. The 
GS due solely to the Dz component is an in-plane 120° 
spin structure invariant under a global rotation around 
the z axis, while a finite Dp component prefers a state 
with a finite uniform out-of-plane spin component, thus 
removing this state from the GS manifold of the Heisen- 
berg Hamiltonian by eliminating the rotation symmetry. 
For Dz > 0, the tilt (/> of each spin from the kagome plane 
caused by Dp is given by tan(2(/)) = 2Dp/{\^J + Dz) 
|48J . which yields (/> = 6° in vesignieite. Interestingly, the 
weak ferromagnetic spin component estimated by NMR 
|27J amounts to 0.05-0.12 /xs, which together with the 
dominant in-plane component being larger than 0.6 ^b 
limits the tilting angle to 3° < < 9°. In the quan- 
tum picture |l0| Neel ordering is induced by Dz only for 
Dz > D^ ~ 0.1. Including a finite Dp in this case leads 
to a weak ferromagnetic moment in the z direction that 
is still linear in Dp/ J as in the classical case [10|, while 
the position of the QCP would be affected by the D^/J 
term as well as linear terms in the AE anisotropy. In 
vesignieite, the condition Dp > Dz, could profoundly af- 
fect the QCP because Dp disfavors spin structures from 
the GS manifold of the isotropic J and should therefore 
be much more efficient in suppressing quantum fluctu- 
ations than Dz- This could explain why magnetic or- 
dering in vesignieite occurs at surprisingly high tempera- 
ture Tn/J = 0.17, despite possessing very similar Dz/ J 
as herbertsmithite. Comprehensive theoretical investiga- 
tions of the general DM-perturbed phase diagram of the 
QKA are thus highly desired. 

In summary, employing the ESR line-width and line- 
shift analyses we have shown that the in-plane DM in- 
teraction is prevailing in the novel QKA vesignieite and 
is most likely responsible for its magnetic ordering be- 
low Tat = 9 K. We have detected intrinsic inhomogeneity 
of the kagome planes, which we attribute to bond dis- 
order, as well as sizable spin-phonon contribution that 
might be related to a lattice instability. Last, we note 
that a preliminary analysis of the ESR line of herbert- 
smithite [5| with the AE model yields anisotropy con- 



stants |A|/J = 0.072 and \E\/J = 0.074 that give an 
"effective" AE anisotropy of |Aav|/^ = 0.06 if averaged 
over the triangle. This value is notably smaller than the 
recent estimate A^^-y/J ~ — 0.1 [30|, which would lead to 
much broader ESR lines since their width scales with the 
square of the anisotropy. Increasing the sensitivity by 
performing single-crystal ESR and applying the above- 
presented analysis is likely the most reliable approach for 
resolving the standing issue of the dominant anisotropy 
in herbertsmithite, which could turn to be the crucial 
milestone in understanding its spin-liquid properties. 
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In an electron spin resonance (ESR) experiment a mag- 
netic system is exposed to the applied static magnetic 
field Bo and an electromagnetic wave with polarization 
of its magnetic field perpendicular to Bq (conventional 
Faraday configuration). Within the linear response the- 
ory the absorption spectrum /(w) is proportional to the 
imaginary part of the dynamical susceptibility []| , 

/•OO 

/(tj)cxx"(q^O,w)cx / dt(S'+(<)5-(0))cxp'"^*/^, 



(3) 
and thus effectively measures spin correlations in the di- 
rection perpendicular to the applied field, {S^{t)S~{0)), 
where ( ) denotes canonical averaging and S" ~ J2i '^f 
is the a-component of the total spin operator. 

Calculation of the time dependence of the spin 
operator in the Heisenberg representation S^{t) = 



ei«*5+e 



i-Ht 



for a general spin Hamiltonian "H is a non- 



trivial problem, therefore, a few approximate solutions 
have been developed. Here, we apply the well-established 
Kubo-Tomita (KT) approach [l| , which relies on dividing 
the spin Hamiltonian into two parts, "H = 'Ho + 'H', where 
the first, dominant part Hq contains only the Zeeman 
term and the Heisenberg isotropic exchange (J) term, 
while the magnetic anisotropy term H' is treated as a per- 
turbation. The latter then determines the shape of the 
ESR spectrum (its position and width), because "Hq pos- 
sessing SU(2) symmetry conserves the total magnetiza- 
tion and therefore leads to a (5-function resonance at the 
field Bo- For finite H' , at high temperatures (T ^ J) the 
KT theory predicts a Lorentzian exchange-narrowed |2| 
ESR absorption line with the full-width-half-maximum 
line width Isll 



(4) 
where 

Af4 = 




[n',s+][s-,n']) 

{S+S-) 
[H - Hz. [W, S+]] [H - Hz, [W, S-]] 
{S+S-) 



(5) 



are the second and the fourth moment of the absorp- 
tion line, respectively, with [ ] denoting a commutator, 
C is a constant of the order of unity (see below), ks 
stands for the Boltzman constant and fiB for the Bohr 
magneton. The expression Q is valid if the magnetic 
anisotropy is small compared to Hq and if spin diffusion 
is negligible, which is generally the case in spin systems 
with dimensionality exceeding one [j|. Strictly speak- 
ing, the ESR absorption line is never truly Lorentzian, 
as all of its moments, given by the spin Hamiltonian, 
are always finite while they diverge for the Lorentzian 
line shape. In systems with strong isotropic exchange 
compared to magnetic anisotropy deviations from the 
Lorentzian shape occur only in far wings of the resonance 
and an approximate line shape that is a product of the 
Lorentzian and a broad Gaussian J3| ex c~^^~-^°^'^^^'' , 
with Be = ks/gt^B \J^iA /M2 being the exchange field, is 
applicable. This then yields C = -\/27r 

Employing this moment approach in the param- 
agnetic limit (T ^ J) when spin correlations on 
neighboring sites vanish, both, the Dzyaloshinsky- 
Moriya (DM) magnetic anisotropy [see Fig. 1(b) of 
the main part] "H^j^ ~ Sfji) ■'-'y ' ^i ^ ^j ^'^"^ ^^^^ 
traceless symmetric anisotropic exchange (AE) in- 
teraction [see Fig. 1(c) of the main part[, "H^ 



E 



(y) 



2A Q? cC 



-f + # ^'5' + 



AE . 

A E\ Qi> ou 

'3 2 / *-^J '-^j 



yield the ESR line width of the same general form 



AB{e) 



kh 



2g{9)fiBJ 



{a + b cos^ 



d cos^ 



(6) 



Here g{0) denotes the 5-factor averaged over the basic 
hexagon of the kagome lattice for a given polar angle 9 
between Bq and the normal to the kagome plane. For the 
DM model, the a, 5, c, d constants are given by the two 
DM components D^ and Dp as [5| 



a - 


= 2Dl + 3Dl, 


b = 


'-2DI-DI, 


c = 


= IQDl + 7%Dl 


d = 


-- IQDl - 2QDl 



(7) 



while for the AE model we derive 



a = 12A2 - 4:EA + 7E^ + {iEA - iA^ + 3E^) cos^6lo, 
b = 12EA - 4A2 + 3E^ + ('12A2 - 12EA - 9E'^) cos^^ 



'0, 



(72EA - 36A2 + 36^^ 



"2^- (36A2 + 36^A + 9E^) cos''6', 







c = 48[68A^ - 36^ A + 33^;'' 

X COS^^o 

d = 48[60-BA - 28A2 + 9E^ + (QGA^ - 24£;A - 36E^) 
X cos^eio - (36A2 + 36SA + 9E^) cos^6io] . (8) 

Here 6^ = 45° corresponds to the angle between the ^ 
axis of a local coordinate system and the normal to the 
kagome plane while the rj axis is parallel to the Cu-Cu 
bond. The direction of the ^ and 77 principal axes of the 
AE anisotropy tensor are determined by the sum of the 
gr-factor tensors of the two sites constituting a particular 
bond [6|. The latter tensors are very close to being uniax- 
ial, with the anisotropy axis pointing towards the apical 



0(1) site [see Fig. 1(a) of the main text], which leads to 
the above-mentioned direction of the ^ and 77 axes. 
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